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NOTATION 


Coefficient 
Projected area of boiler rear wall, sq. ft. 
Area of burners = nvd*, sqe ft. 

Exponent of Reynolds Number 

Constant of integration 

Thickness of boundary layer 

Constant 

Specific heat, BTU/lb. - degrees F, 

Furnace depth - ft, 

Burner or orifice diameter - ft. 

Mass flow in ‘units of 1b/sq.ft. - hr. 

Heat transfer coefficient, BTU/hr. - degrees F. = sq. ft. 
Thermal conductivity, BTU/ft.- degrees F, - hr. 
Subscript denoting mean value 

Number of burners 

Two dimensional quantity of flow 

Heat transfer per unit time, BTU/hr. 

Air to fuel ratio by weight 

Distance from side wall to burner 

Effective cold heat absorption surface 
Projected heat transfer surface in rear wall, sq. ft. 
Absolute temperature, degrees R, 


Absolute temperture of cold surface 


' Furnace exit temperature from radiation calculation 


Effective flame temperature 


Film temperature 
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Steam or boiling water temperature 
Actual exit temperature from furnace 


Coefficient of convection heat transfer, 
BTU/hr. = degrees F, = SQe tt. 


Velocity in x direction 

Velocity in y direction 

Free stream velocity along rear wall 
Weight of air, lbs/hr. 

Weight of fuel, lbs/hr. 

Weight of gas, lbs/hr. 

Complex potential function 

Complex plane, x+i1y 


Nusselt Number = hD 
ae 


Reynolds Number = GD = Le '2 
ZE ZO 


Prandtl Number = Cy x 
k 


Temperature difference between free stream and wall 
Log mean temperature difference 

Constant for heat transfer equations 

Dynamic viscosity 


Absolute viscosity 
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INTRODUC TION 


Since the end of the nineteenth century, the marine 
boiler, both merchant and naval, has gone through an extensive 
period of evolution. The firetube boilers,such as the Scotch 
boiler,have given way to watertube boilers of either the header 
or drum type. In more recent years, the integral furnace two- 
drum boiler with integral uncontrolled superheater has gained 
considerable favor due to its simplicity, relatively light 
weight and low cost. In order for these boilers to meet the 
specified superheated steam temperatures without control, it 
is necessary to pay considerable attention to the proper design 


of all heat transfer surfaces in the boiler. 


All boilers possess some sort of combustion chamber or 
furnace for burning the fuel. The gaseous products of combustion 
pass over the generating and superheating surfaces prior to 
entering any other heat recovery equipment and exhausting to the 
atmosphere. It is obvious that in order to effectively evaluate 
the heat transferred in the generating and superheating sections 
and thus the superheat temperature, the temperature of the gasses 


leaving the furnace must be accurately determined. 


The method of calculating the furnace exit temperature is 
one of equating the heat given up by the gasses to the heat 


absorbed by the furnace surfaces and solving for the exit tem- 


(1), 


perature Numerous descriptions covering the several methods 


of attack on this problem may be found in the current litera- 


ture !192s32ls556) | In general, the heat transferred to the 
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furnace surfaces by convection has been neglected and the 
problem has become one of radiation only. Until recently 
good results have been obtained using these methods, whether 


empirical, theoretical, or a combination of both. 


In recent years the severe space and weight limitations 
imposed on naval boilers in particular, combined with increased 
demands on the boiler, have forced the designer to use heat 
release rates in.excess of 500,000 BTU/ou. ft. In order to cool 
the surfaces of the furnace, water-cooled walls, floors, and 
overheads have been adopted. These inovations were adequately 
handled as cold surface when using the conventional method of 
solution of the furnace problem. However, as increased space 
limitations have dictated decreasing the depth of the boiler, 

a discrepancy arose between the calculated and measured super- 
heater outlet temperature. The error can be traced to too high 
an estimate of furnace exit temperature. It appeared that heat 
was being absorbed in the furnace by some means other than 
radiation. Examination of the designs involved showed that they 
were fitted with shallow furnaces and watertube rear walls. 

The additional means of transfer was by convection through 


impingement of the gaseous flame on this rear wall. 


Since no data on this means of heat transfer was available 
in the literature, early evaluations were by "rule of thumb" 
empirical means. Bethge and mieeda’ carried out an in- 
vestigation on one particular boiler. Holmboe and Hove ( 8) 


investigated impingement on a flat plate with the use of 
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flat plates as the surface instead of tubes. It is because 


of this lack of adequate basic data that the present investi- 


gation was undertaken. 
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GENERAL CONSIDERATIONS 


The problem of heat transfer by fluid impingement on a 
boiler rear wall resolves into a problem of forced convection. 
Since the heat is conducted through a turbulent boundary layer 
or film, the thermal conductivity of the fluid k will be a 
factor. ‘The film thickness depends on the mass velocity through 
the burners G, the distance between the rear wall and the 
burners or furnace depth D, and the viscosity of the fluid 2 5 
For a given quantity of heat transferred, the temperature of the 
fluid stream depends on the specific heat Cp. Letting ¢ repre- 
sent a function, the following equation for the coefficient of 
heat transfer h may be written: 

Vicaid GID 1G. Qa eas) ace 0 “« @ cane 
This may be written in an infinite series, all terms of which 
will have the same form. Therefore, taking the first term only: 

ee Grew ok ass ee hare a 
where the factor o_is an arbitrary constant dependent on the 
geometry of the furnace and the boundary conditions. Substi- 
tuting the dimensions and solving simultaneously for the ex- 
ponents in terms of a and e gives the following equation: 


LD = (Ss (ee. ame a ee 


In more familiar terms: 


On 
Ny =o Nee Ne avis. enter le’ Pewee en) 


Ny 
The problem now remains to find the values of the exponents 
@ and 6 and the constant @. . Examination of the geometry of 


the fluid flow in the furnace reveals that the fluid leaves 
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the burners ina cone shaped stream directed at the rear wall. 
However, in order for continuity to be satisfied, the fluid must 
eventually flow parallel to the rear wall and out through the 
tube bank. Therefore, the convection heat transfer to the 

rear wall is conducted through a boundary layer of fluid stream- 
ing parallel to the surface. For this reason, the constant © 
must contain those parameters that have a significant effect 


on the thickness of this layer as well as its velocity. 


For a given weight of gas flowing, G will be dependent 

on the size and number of burner orifices. Instead of using 

the number of burner orifices, the ratio of projected rear wall 
area to burner area Ap/Ao will be more flexible while still 
defining the significant factors in the geometry. Also evident 
is the fact that increasing D will increase Nre proportionately; 
however, the heat transfer will obviously not increase due to 

the decrease in fluid velocity as depth increases, Therefore, 
the ratio of furnace depth to burner orifice ene ets was 
chosen as a significant parameter. This ratio, together with 
Ap/Ao, fix the significant features of the geometry. Therefore, 


the constant o will be a function of these parameters D_ and aa 
a Oe 


Holmboe and Hove (8) conclusively demonstrated the validity 
of this approach. Using two similar models having a linear 
size ratio of 2:1, they were able to measure the same value of 
Ny for equal values of the dimensionless ratios GD Se 
Therefore, they justly concluded that extrapolation of model 
results to full size by using the dimensionless Bquation was 


a valid approach. 
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METHOD OF APPROACH 


Since the geometry of the problem 4s so complex and the 
variables so numerous, the easiest and by far the most depend- 


able method of finding a, 6, and Q was by experimentation. 


The method used consisted basically of testing a model 
boiler on which the parameters Ap and D could be varied. 
The actual phenom€num of heat ee was the reverse of that 
found in an actual boiler in that the rear wall heated the 
gasses flowing through the furnace instead of cooling them. 
Air was blown through the burner orifices of the model, im- 
pinged on the rear wall, and then passed out the stack. By 
measuring the weight of air flowing Wa, the inlet and outlet 
temperatures Ty and To, the heat transferred could be computed. 
From this, Nusselt Number and Reynolds Number were calculated 
and plotted for each change in geometry. The apparatus was 
tested over a range of 5 values of_Bfrom 3 to 10.5 and over a 
range of }} values of Ap_ from 11.47 to 45.8. The weight of air 
flowing Wp varied (ken atous 260 to 1200 lbs/hour. In all,data 
at 567 points was taken. For a more detailed description of 


the apparatus and testing procedures, see Appendices A and B. 


In the early stages of testing, it was necessary to run 
many of the tests twice or even three times to be sure that 
the data was reliable. However, as the testing procedure be- 
came more refined, this was no longer necessary. By insuring 
that steady state conditions existed in the model, reliable 


and repeatable results were obtained in almost all cases. 
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Figures 1 through 5 are plots of Ny, vs. Neo for the 
five values of uw and air flowing through plain orifices. 
Figures 6 through 10 are plots of Ny, vs. Neg for the same 


five values of 3 but with the air flowing through model burners. 


In addition to the experimental approach an attempt was 
made to develop a theoretical analysis of the problem. This 
analysis consisted of senaideeine the two dimentional problem 
only. Using conformal mapping techniques, the fluid flow in 
the furnace was mapped giving the free stream velocity distri- 
bution across the rear wall. Then considering the problem as 
one of heat transfer through a turbulent boundary layer on a 
flat plate, a relationship for Nusselt Number in terms of the 
Reynolds Number using the velocities previously calculated 
was derived and the results of this derivation compared with 
the test results. The differences in results are covered 


in the discussion. 
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ANALYSIS OF TEST RESULTS 


Figures 1 through 10 are plots of Nyy VSe Neg plotted on 
log-log paper. It is evident that for each value of 5 and 
pe» the test spots plot in straight lines. Therefore, Nyy is 
definitely a function of Nag to some power a3 and a will be 


equal to the slope of the lines. 


With the equipment available for these tests, it was im- 
possible to vary the Prandtl Number of the gas. Therefore, 
the value of the exponent 4@ is assumed to be 1/3. This is 
in agreement with the equations for forced convection heat 
transfer over various shaped surfaces as given by McAdams(3)} 
and Jakob '9»10,11) The differences in most of these equations 


lies in the exponent of Nre and in the constant. 


From Figures 1 through 5 for the case of air discharging 
through plain circular orifices, the value of a is 0.87 for 
all curves. There are some slight variations in slope, but 
they are small enough to be considered negligible. Therefore, 
for fluid discharging through circular orifices and impinging. 
on the boiler rear wall, the equation of convection heat 
transfer is: 37 

ha A (SP) (ae) ee ne ae 
Where the values of A, are given in Table 1l, and Cp, <, and 


k are evaluated at the film temperature. 


For the case of air discharging through burners, analysis 


of Figures 6 through 10 yields the value of a as 0.90 for B a3), 
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5, 7, and 9. However, at 2 210.5, the value of a increases 
to 1.0. Therefore, the following are the equations for con- 
vection heat transfer with fluid discharging through burners 


and impinging on the boiler rear wall: 


EO. aN 

for os 9 aL = A eee) (ee) © Po aie eer a) (6) 
: HANS. 
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where the values of ~, are given in Table 2 and Cp, iy and 


k are again evaluated at the film temperature. 


Table 1--- Values of o( for Equation 5 


Air Discharging Through Plain Circular Orifices. 
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Table 2--- Values of for Equations 6 and 7 
Air Discharging Through Typical Marine Burners. 
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THEORETICAL ANALYSTS 


Considering the two dimensional case only, the plan 
view of the boiler furnace may be considered as a semi- 
infinite strip in the z plane ‘12), The burner is assumed 
to be a source in the front wall a distance r from the side 
wall. The effective width of rear wall is assumed to be 2r. 
Thus, the boiler furnace with one central burner orifice can 


be represented by the following figure: 


Zz CRANE 


Using the Schwarz-Christoffel heaven! 


» the solution 
for the free stream velocity along the rear wall may be obtained 


in the following manner: 


The flow in the t plane may be represented thus: 
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where H is at t= = <1, Gat t = 1, and Fatt= k 
and k > l. 

The deflection angles at H and G are /2 and at F it 
is zero. Therefore, the Schwarz-Christoffel theorem 


may be written as follows: 
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Applying the boundary conditions as follows: 
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Since the total strength of the source is -Q, in the 
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Therefore, <= © © is a stagnation point also. 


To get free stream velocity along rear wall, 
Avo 


Where 7, and wW 


4 are velocity components in the x and 
y directions respectively. 
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Therefore, equating real and imaginary parts, along 


the rear wall: Au = S 
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It will now be necessary to arrive at a suitable equation 
for N, using this velocity. Due to the complexities of the 
geometry, attempts to analyze the problem in three dimensions, 


or even allow for the fact that the surface is composed of 


tubes, proved fruitless, Therefore, it became necessary to 


simplify the problem. 


As an approximation one can consider the boiler rear 


wall as a flat plate with a turbulent boundary layer,’ 
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The origin of the axis system chosen will be at the 
stagnation point in the corner between the side ano rear 
wall. To simplify the problem greatly, only the two dimensional 
case will be considered. It can be seen that for unit width 


in the z direction: 


Es dg dy where Ag. is heat trans- 


ferred to the wall over a unit length in the x direction. 


> lo 
ee = a dq ay = ig a 72 Ce SAP du. 
where \ is the velocity of the fluid at y 


Avis the temperature of the fluid above the wall 


temperature. 
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Combining terms and dividing by ax gives: 
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However, von Karman ‘9 ) gives the velocity distribution as: 
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N= Ng (+) ? 
and Rube sin‘10) » Simplifying the equation for b by 


omiting the buffer layer and assuming a somewhat thicker 
boundary layer, gives the following relationship for b: 
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where the subsceipt x denotes the quantity evaluated at x. 
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Selecting values of depth between r and lr, and substi- 
tuting values of x from 0 to 2r into Equation 8 gives values 
of |4,x| . Raising each value of [Yo*4 to the .8 power 

¥ “We 
] . = 
and taking a Simpson's mean gives Yo, x] ae om f=" ri 


where Cy is given below: 
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Noting that (221 > oar and substituting the above 


values of Wo, a into Equation 9 gives: 
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Since r =D, D/2, D/3, D/k respectively, the following 


equation may be written: 
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where , is given by Table 3. 
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Table - Values of Xx for Equation 10. 
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DISCUSSION OF RESULTS 


In comparing the results obtained with circular orifices, 
burners for D< 9,and burners with 5 > 9,care must be taken 
to note that each case has a different power of Nre as well 
as different values of 4 . For example, using air with 
Nro = 300,000 and N,, ® -7 gives: 


Circular orifice —~— Re aa 
Burners D = 9 —_- N Bars -& eGha ; BOO) 
. a Nu 
> PD = 6 (2,66,O00C 
Burners D > 9 Nyy oa ‘ ) 


Although the values of Nyy will diverge somewhat due to 
changes in Np, » the values for the above case are plotted in 


Figure 11 and several general trends are evident. 


Nyy with circular orifices is roughly 50% greater than 
with burners. Since the use of burners imparts a whirling 
motion to the fluid, the sideward velocity is considerably in- 
creased while the forward velocity is decreased. Therefore, 
the convection of heat by efflux into the boundary layer will 


be proportionately reduced and the burners will give a lower 


value of heat transfer. 


Another interesting aspect 4s that beyond 5 of 9 the 
curves of Nyy for burners falls off sharply. It is apparent 
that beyond this point the free area through the screen tubes 
is large and the forward velocity small so that the fluid 


fails to reach the rear wall in sufficient quantity to prevent 
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some stagnation from occuring. The net atfect is a resultant 
decrease in heat transfer. With circular orifices, however, 

the major component of velocity is forward; and the fluid reached 
the rear wall in all cases within the ranges tested. From this 
dropping off, it may be assumed that in actual boiler furnaces 
the convection heat transfer to the rear wall is negligible 
beyond 5 *» 9. This is substantiated with observations on 


marine boilers where p of 9 would mean a depth of about 12 feet. 


Although the curves for = of 22.9 «045.8 cross the other 
curves for the orifice case, the general trend of both the 
orifice curves and burner curves is the same. Also, With 
reference to Figure 12°'Ft is evident that as AB decreases, & 
first increases to a maximum value and then decreases again. 
For the circular orifices the value of Ap for maximum a. 
becomes larger as depth increases, reat be interpreted 
to mean that at greater depths maximum heat transfer will 
occur with smaller orifices. For the burners, all the curves 
are similar and maximum <q, is at = of about 13. Beyond that 
point there is no advantage in adding more burner area. 

A search of the current literature revealed that 
McAdams (3) gives the following equation attributed to Colburn 
for heat transfer through a turbulent boundary layer with flow 


parallel to a plane surface: 
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This equation is of the same form as Equation 9 except 
that the Prandtl Number, [Ces i, | is to the 
1/3 power. In addition, McAdams gives the following relation- 
ship attributed to dJurges for flow of air at room temperature 


parallel to a vertical copper plate: 
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where aia) is expressed in “°°/@- FT *- Dec F. wamac AT. 


and V* 4s in ft/sec. The factors for Equation 12 are given 


in Table l. 
Table - Factors for Equation 12 
SE een ||| Vi <te evsec| 16 race < vid 00 evsee 
SORT oie sures 
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Tt is clear that the value of 4, is dependent on (Vv) 


1 with vt < 16 ft/sec. 


when wes 16 ft/sec. but 1s dependent on V 

With the burners at = = 10.5, Equation 7 shows the same 

dependence on the velocity, rather than the velocity to the 

-9 power. AS previously stated, some stagnation must occur on 
b 


the rear wall beyond 3 - 9, and thus the mean velocity would 


be very low - even below 16 ft/sec; 


' Since very few tests have been run on heat transfer by 


fluid impingement, there is very little data available for 
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comparison purposes. Friedman and Mueller! 4 ran some tests 
on air impinging on a horizontal heated plate and obtained 
the following equation: 

dr ec aed 
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where C, was dependent on the plate and the spacing of the 
discharge orifice from the miller C, varied from .011 to 
e321 depending on the geometry. The geometry of these tests was . 
considerably different from that of a boiler furnace so no 


comparison can be made, 


(8) 


Holmboe and Hove give the following equation for 
_ impingement on a flat copper plate in the same boiler model 


used in these tests: 
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The value of of -0258 is slightly higher than the values 
of .0095 to .021 found in the present tests for air discharging 
through orifices. This discrepancy is due to the lower ex- 
ponent of No found by Holmboe and Hove. The lower value of 

the exponent is believed due to inability of the equipment to 
Bisintain steady state at high rates. This caused their plots 
to curve over, and thus the resulting line had a lower slope. 
This equation gives a walue of Ny, of 738 for the case of 
Non = of and Nee = 300,000. This is approximately the mean 
value of the curves of Nyy for the orifice oase as plotted in 
Figure 11. It is of interest that the use of tubes instead of 
a flat plate made very little change on the heat transferred. 
This would support the use of the projected rear wall heat 


transfer surface in the calculation of h in Appendix E. 
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In analyzing the theoretical solution, no direct com- 
parison can be made since Np, is actually based on the 
velocity along the rear wall rather than the mass flow through 
the burners. Since G = Wg and A, is zero for a source, 

G has no meaning in this BeAceuen as it would be infinite. 


However, since the Reynolds Number is represented by (oa 


AL 


it is of interest to note that oS is the mass flow along the 


aes 
rear wall. Geese” equaled == where G is mass flow 

through the biuners, a direct comparison could be made. 

Since for a constant value of Reynolds Number, Nyy is pro-=- 
portional to %, , aA, is plotted in Figure 11 for the sake 

of examining the trend of the curve. In the model, r was 6" 

and d was 2"3 therefore, for comparison purposes D may be assumed 
to equal 3 De. It may be seen that the shape of oe curve is 
similar to we burner curves with the exception that the drop 

off beyond D = 9 is not included in the theory. Tt is inter- 
esting to aii that the curve is steep at shallow depths and 


levels off as depth increases, 
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APPLICATION TO DESIGN 


As previously stated, the method of calculating furnace 
exit temperature is one of equating the heat given up by the 
gasses to the heat absorbed by the furnace surfaces and solving 
for the exit gas temperature. The rigorous solution of the 
problem would involve equating radiation plus convection to 


the heat given up by the gas using the following equation(2), 


4 = 
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where S, is cold surface, FF, is the combined shape emissivity 

factor, Tp is effective flame temperature, T,. gas side tem- 

perature of the cold surface, U, is the coefficient of convection 

heat transfer, T,, is the mean temperature difference, Sw is 

the projected heat transfer surface of the rear wall, Wp is 
weight of fuel, R is air to fuel ratio by weight, LHV is the 

low heating value of the fuel, Ge is the enthalpy of the fuel 

above to, (tg - to) is the temperature of the air above to» 

Cy is the specific heat of the incoming air, and an, is the 


sensible heat of the flue gas at the furnace exit temperature Tre 


However, the question remains as to what value of AT, 
should be used in the convection problem. The flame temperature 
Tp is too high a gas temperature since the gasses in the vicinity 
of the boiler rear wall have been found to be considerably 


cooler(15» 16) In obtaining the data presented herein, the 
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gas was assumed to enter the boiler front at some constant 
{nitial temperature. In the actual furnace, combustion takes 
place after the fuel and air enter the furnaces; while in the 
test apparatus, combustion was assumed to have taken place 
prior to blowing the gasses into the furnace. Therefore, in 
order to use the data as presented, the assumption must be 
made that combustion has taken place and the gasses have 
assumed some initial temperature prior to impinging on the 


rear wall. 


The easiest way to handle this is to first calculate 
the furnace exit temperature Tp using only the radiation 
portion of Equation 15. The gasses are then assumed to impinge 
on the rear wall at this temperature Tre Using the following 


equation, a corrected temperature of the gas leaving the 


furnace T, may be computed: ‘ 
Us Su, 
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where T, was previously calculated.and re is the outside 
temperature of the rear wall tubes,and may be considered 
equal to the temperature of the boiling water within the tubes 
Tg, Since the coefficients of heat transfer through the metal 
and to a boiling liquid are very high. Sw is the projected 
area of the rear wall heat transfer surface, We = Wa+ Wes 


Cy is the specific heat of the gas evaluated at the average 
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‘ 5 and U, is given by the following: 


gas temperature 


alle. os es +~ ae ae \ 


Oe a Wea tee, | 
For change of state, the coefficient of heat transfer through 
the water film inside the tubes hy is very large; so 1 may 
be neglected. Likewise, k, the thermal conductivity of the 
tube wall, may also be neglected unless the wall thickness is 


large. Therefore, U. ~% hh » and h may be found from 


Equation 6 repeated below: 
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where Cy » “~~ , and k are evaluated at the film temperature 


Tp and Tr = Tg+4Tm wherefTm is given by the following 
ad 


equation: 
Ge = 
ISS W = ee e« © © @ (18) 
Me (ee Res 
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Values of O. are plotted in Figure 13 vs. D/d and Ap/Age 
In computing G, the vanes in the registers are not allowed 


for; therefore: 
Wa* 
a ee eo ae we o CSD 


where Wy is the total weight of gas and A, is the total 


burner area. 


By estimating T), the film temperature and average gas 
temperature are computed. Using these, U, and Cy are determined 
and T,; solved for using Equation 16. With this new T,, the 


process is repeated solving for another value of T. Usually 
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accuracy for engineering calculations. 


This correction for rear wall convection may be considered 


ligible with D/d greater than 9. 
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CONCLUSIONS 


One of the basic assumptions made in analyzing this 
problem is that the gas enters the boiler front at some 
{initial cori@ition. In actual practice, oil and air enter 
through the atomizers and air register respectively; the 
oil ignites, and the burning gasses expand outward. The 
total weight of gas is equal to the weight of air and oil to- 
gether. This doesn't change in the combustion process; but 


the velocities do. 


In spite of this, the data presented for the fluid dis- 
charging through burners should give a good approximation of 
the heat transfer by convection in the boiler rear walls. 


The method of applying this data is given in tue pr&teeniwe secrTiow. 


In general, the following conclusions may be drawn in 


regards to the convection heat transfer in boiler rear walls: 
1. At depths beyond D/d = 9 the heat transfer is negligible. 


2, The heat transfer is a function of D/d and Ap/Ao, and 


the use of these parameters in testing was correct. 


3. There is no increase in heat transfer by reducing 


Ap/Ao beyond a value of about 13. 


lh. In computing the heat transfer, the pro jected area 


of the heat transfer surface may be used. 
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RECOMMENDATIONS 


The following recommendations may be made as a result 


of this study: 


1. Tests should be run on models with the tubes at 
various spacing and backed by both tangent and "T" tile 
to determine if the use of the projected area would be 


justified in these cases, 


2. Further small scale testing such as this should 
be done using liquids as the fluid. This would allow more 
accurate measurement of velocity, weight of fluid, and 
temperature since radiation would be virtually eliminated. 
The use of liquids would also allow the Prandtl Number of 
the fluid to be easily varied. 
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APPENDTX A 


DESCRIPTION OF APPARATUS 


In order to effectively evaluate the convection heat trans- 
fer on the boiler rear wall, it was necessary to simulate with 
models the fluid flow and heat transfer phenomina that conus in 
a boiler furnace. To accomplish this, a model boiler was used. 
Air supplied by a turboblower through a flow meter was blown 
through the model. The heat transfer was accomplished by re- 
versing the normal boiler process and allowing the heated rear 
wall to heat the air. By measuring the weight of air flowing 
and inlet and outlet temperatures, the actual heat transferred 


to the flhid could be computed. 


Therefore, with reference to Figure A-l, the testing 
apparatus consisted of a model boiler, steam generator, turbo- 


blower, flow meter, and associated valving and instruments. 


The model boiler (see Figure A-2) was a modification of 
the 12 inch model used by Holmboe and Hove | & ), It consisted 
of a plywood box similar to a boiler furnace with an adjustable 
front, tubular rear wall, and a row of dowels representing the 
screen tubes in a marine boiler. The furnace front was a@ mov- 
able copper plate containing the burner orifices. The scale 
chosen was 1/8 inch = 1 inchs; so that the two inch orifices 
represented 16 inch diameter burners. By moving this plate, 
the furnace depth could te varied from 6 to 21 inches. Four 
different orifice plates were constructed with 1, 2, 3, and 4 


orifices respectively. The orifices were arranged so that the 
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center of area of the orifices coincided with the center of 
the plate. In addition they were cut 32" between centers to 


correspond with normal boiler practice of 26 inches. 


The tubular rear wall (see Figure A-3) was a series of 
thirty-two 3/8" 0.D. copper tubes brazed into an upper and 
lower header. The lower ends of the tubes were curved through 
Slightly less than a right angle to give flexability for thermal 
expansion. Sufficient length was allowed so that a 12 inch by 
12 inch flat surface of tangent-=tubes was exposed to the im-=- 
pinging fluid. The upper header was fitted with a steam inlet 
and a combination vent and pressure tap. The steam wes dis- 
tributed in the header through a length of 3/8" 0.D. copper 
tubing, plugged at the end, and deitied with 1/16" holes whose 
combined area was equal to the cross sectional area of the inlet 
line. The lower header was fitted with a condensate drain. 
As installed in the model, the tubes were backed by 1 inch of 
asbestos cement, and the headers insulated with rock wool pad- 


ding. 


Steam for heating the model rear wall and to run the 
turboblower was supplied by a Clayton Steam generator. Steam 
to the rear wall passed through a needle valve into a moisture 
separator. A small tangential hole admitted the steam with a 
swirling motion, separating the entrained water which drained 
out the bottom. The dry steam passed out the tap to the upper 
header and into the tubular rear wall. The condensed steam 


drained into the lower header and out the condensate drain. 
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+ The turboblower was an Allis-Chalmers aircraft turbo- 
supercharger with the nozzle blocks modified to take steam. 
Saturated steam at 110 psig, as provided by the steam generator, 
passed through a regulating valve and a throttle valve provid- 
ing a nozzle chest pressure of 50 psig. A bypass was used to 
direct any excess air back through the turbine and out the 


exhaust. 


The flow meter was a concentric, square-edge orifice meter 

fitted with flange taps. The meter was installed in sufficient 

straight length of 3 inch standard pipe to insure accurate read- 
ings. A Deni tach diameter orifice was used, and the pressure 

drop across this orifice was read directly in. inches of water 


on a differential manometer. 
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In addition to the above, manometers were installed to 
_ measure positive steam pressuré in the model rear wall and air 


pressure in the inlet chamber of the model boiler. 


Two mercury thermometers were installed, one in the air 
inlet chamber and one in the stack, to measure the AT of the 
air through the boiler. They had a range of 51° Cc. and were 
graduated in Gal Cx allowing the temperature to be read quite 


fo) 
accurately to the nearest .02 C. 


A small balance and a timer were originally used for 
Peesuring and timing the steam condensed in the rear wall. 
However, this was eliminated when good agreement could not be 
obtained between the heat of condensation and the heat absorbed 


by the air. 
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Model burners similar in geometry to those used on marine 
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boilers were constructed. They could be attached to the ori- | 
fice plate in the furnace front and imparted a swirling motion _ 
to the air as it entered the furnace. Details of the burners 


are shown in Figure A-l. 


A description of the test procedure and method of oper- 


ation of t he equipment is given in Appendix B. 
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Figure A-5 View of Model Boiler Taken from Rear Wall 
End 
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Figure A-7 View of Model Boiler With Doors Open 


Figure A-8 Burner Orifice Plate with Four Burners 
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Figure A-9 Steam Driven Turbosupercharger as 
in the Tests 
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APPENDIX B 


TEST PROCEDURE 


The procedure for preparing the model and running a test 


is as follows: 


1. While the steam generator was coming up to pressure, 
the model was made ready for the test. The desired orifice 
plate and partition necessary for a given depth were selected 
and assembled in the model with the orifice plate against the 
screws in the floor. After insuring that the split rubber 
tubing on the partition and the rubber tubing seals around 
the doors and orifice plate assembly were properly located, the 
doors were closed and tapped with a hammer to set the points on 
the orifice plate assembly into the wood. The trunk snaps were 
then closed, squeezing the rubber tubing and making an air tight 


seal. 


2. At this point all the manometers were adjusted to zero 


and the two thermometers checked for identical readings. 


3. When the steam generator was up to pressure, the 
turboblower was warmed up and the speed slowly increased to 
maximum speed. The bypass control, in the combustion chamber 
on the turbine, was closed and the air regulating valve opened 
Wide (see Figure A-1). Steam was next cut into the rear wall; 
and the vent valve, water trap drain valve, and condensate drain 
valve adjusted to allow a wisp of steam to flow out. This in- 
sured that the rear wall was completely filled with steam and 


no condensate or air were trapped in it. The steam pressure 
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in the rear wall was adjusted with the needle valve to a posi- 
tive pressure of one inch of mercury. This pressure was main- 
tained throughout the testing. The apparatus was allowed to 
run at this maximum rate for about 20 to 30 minutes until it 


was completely warmed up. 


h. During the above warm-up time, the barometric pressure 
was recorded as was the wet and dry bulb temperature which was 
taken with a sling psychrome ter. These readings were also taken 
and recorded at the end of the test. The average were used to 
determine the specific humidity using the psychrometric charts 


: ay 
in Ellenwood and Mackey! ate 
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5. After the apparatus was warmed up, the bypass was 
opened to the turbine and the weight flow of air adjusted by 
closing the air regulating valve. Originally the weight flow 
was adjusted by throttling the steam to the turbine, but this 
proved unsatisfactory as the blower tended to hunt. Therefore, 
in order to maintain a steady state, the turbine was run at top 
speed at all times and the weight flow of air regulated by by- 
passing the excess. This method of control proved quite satis- 
factory. , 

The steam préssure in the rear wall was then adjusted and 
the transfer of heat allowed to reach a steady state. This 
could be determined by observing when the inlet and outlet | 


temperatures remained fairly constant. 


6. When steady state was reached, the burner orifice and 


flow meter heads were recorded in inches of water. The first 
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temperature reading was then recorded, followed by the second 
and third at two minute intervals. The flow meter head was 
again checked, and if unchanged, the next weight flow would be 
tested, If a change in flow was noted, the run would be re- 


peated. 


7. After adjusting to the next weight flow, and while 
awaiting steady state, the average inlet and outlet temperature 
was calculated. Using this T, and To, 4 T through the model 
was recorded. The value of AT WARS was then plotted 
against AH - This gave a good Shee ten the consistency 


of the data from point to point since: 


Wee ve VAS 
ae a Q 
Bus Qe We Co AT and Sa = Se 
Therefore , Dinah AA ASN ae OBE 
Ae ATan 
or since Ce is nearly constant, yl Asean ae 
[SSIS 3 
r and since Wereen ee , Nog ~ Ay are 
t 
So Se ANS) 
Also Nee = wez= and G — 7a 
Therefore | Nie ~ af OWE 


Since Nyy, was assumed to be proportional to NRe to some 


power, then on a log-log plot ISOS EN SV vee =/ Ade 
should plot as a straight line. This proved to be an excellent 
way of checking for points in error while still testing. A 


sample plot is shown in Figure B-l. 
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The following pages contain the data as recorded in 


the actual tests: 
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APPENDIX _D 


‘ 


CALCULATION OF WEIGHT FLOW, SPECIFIC HEAT, 
VISCOSTTY, AND THERMAL CONDUCTIVITY. 


The weight of air flowing is dependent on the following: 


1. Geometry of the fluid meter. 

2. Head loss across the fluid meter. 
3. Specific humidity of the air. 

lh. Air temperature. 


5S. Air pressure. 


Factors 1 and 2 are considered when using the method of 
calculating fluid velocity, as recommended by The American 
Society of Mechanical icipinseras — ). Mass flow, however, is 


also dependent on the specific volume of the fluid. 


Using the gas equation, {.@.s, PW = RT, it is evident that 
v is dependent on R, T, and Sb. The specific humidity of the 
air has the effect of changing the molecular weight of the 
fluid and thus the gas constant, R. Since the value of R is 
not very sensitive to small changes in specific humidity, an 
average value of 50 grains per pound of dry air was assumed 


only for the purpose of calculating the weight of fluid flowing . 


Since the air temperature varies considerably from day 
to day over the year, the temperature of the air had to be 
accounted for. An assumed temperature of 100° Pe wes used in 
calculating the weight flow curve, and corrected to the actual 


temperature by the following relationship: 
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Since A4.~ and c~Vau, then ¢~/v but G-= & and v~r 
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From some early testing it became evident that the up- 
stream pressure varied from run to run. The maximum pressure 
was less than 30 inches of water or about 7% of atmospheric 
pressure. Since the velocity of flow c was proportional to 

p, the maximum error introduced by assuming atmospheric pres-~ 
sure is about 34. ‘However, when metering an expansible fluid, 
an adiabatic expansion factor Y must be introduced “8) , This 
factor, less than unity, will reduce the flow. Considering the 
maximum pressure of 30 inches of water, the pressure ratio, | 
r=P2 = .93. | 

PI e 


From Figure 97 of Reference 18, using <= -689, the adi- 


abatic expansion factor Y, = .975. 


Therefore, by considering the pressure at atmospheric, 
1.e., r - 1.0, Y; = 1, the maximum error is: 


1.035 (.975) = 1 = .009 or 9/10 of 1%. 


In addition, since both Nusselt Number and Reynolds Number 
are directly proportional to Wg, both values have the same 
maximum error or .9%. Since Ny, 1s proportional to Npeo, 
the maximum error introduced in the actual plots of Ny,, vs. 

Nhe is of the order of (1.009) °87 - (1.009), or .12%. Since 
this is the maximum possible error and decreases with decreased 
weight of fluid flowing, it is negligible. Therefore, the 
pressure was assumed at 14.7 psia and no correctien made for 


changes in pressureé.. 
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Values of W, corresponding toAHp values of 0 to 30 
inches of water were calculated and are plotted in Figure Dee. 
The following is a sample calculation showing the method used: 
C = Cava 
where C is Seaéed ty of fluid flowing in ft/sec, is head 
loss through the flow meter in feet of fluid flowing, and 
C4 is the coefficient of discharge. Since A Hp was measured 
in inches of water, then: 
ey ee «= x NN Se: Oct 
where \Y is the specific volume of fluid flowing in cubic 


feet per pound. 


Therefore}; 


C= Cy fa aaa «ass ed eax VOTE 
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Assuming an average specific humidity of 50 gr/lb. dry 
air, a temperature of 100° F, and a standard pressure of 


14.7 psia, and using: 


M, = molecular weight of air = 28.970(19), 


a 
MH>o = molecular weight of water vapor = 18.016! 19), 
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Again at the conditions assumed: 
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Where 2 is D2/D1>5 D, is inside pipe diameter, and Do is 
orific diameter, and K is the coefficient from Table 12,ASME‘18) , 


ot RES 2 es BSG 
‘ eX 3.00% 


Ch= WK -/ Gew)t= Sst 


Values of K interpolated from Reference 18 are plotted in 
Figure D-l. 


Assuming OH, = 16 inches, then /AHp = 4.0. Estimating 


C4 as .609 gives the following: 
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The values of specific heat were computed for a pound of 


wet air by the following equation and are plotted in Figure D-3: 
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where Cee is specific heat of water vapor (19) , 
Ce, 18 specific heat of air(19), ana 
x is specific humidity in gr/lb. of dry air. 


Viscosity and thermal conductivity are taken directly 


from the table for dry air{19) and are plotted in Figure D-l. 
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APPENDIX E 
CALCULATION OF NUSSELT NUMBER 
an 
REYNOLDS NUMBER 


For each test spot $he Nusselt Number and Reynolds 
Numbér was calculated and plotted in Figures 1 through 10. 
‘The heat transferred was equal to the weight of gas times 
the specific heat times the change in temperature. Using A Hp 
awl Tlave from the test data, W, was determined from Figure D-e. 
Also usingAT, the values of film temperature Tr and the mean 
temperature difference AT, were calculated using the following 


relations: 


We = oR +87 AWD Aaa ies *S 
4, ee 


The heat transfer coefficient h equals the total heat 
transferred divided by the area and the mean temperature 
difference. The projected rear wall area Ay was used as the 
heat transfer surface giving the following equation for h: 


where the specific heat Cy was evaluated at the bulk temperature 
of the gas T, and corrected for moisture content. Values of 


Cp were taken from Figure D-3. 


From Figure D-, values of k and < at the film temperature 


were taken. Using Nyy = bp the value of Nusselt Number was 
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readily calculated. 


Reynolds. Number was calculated using the following 


equation: 
Ne. = GYD 
e AL 
where G = W, and A, was the total burner area and W, the 


Ko 
total weight of air flowing. 


The actual calculation of Nusselt Number vs. ‘Reynolds 
Number was done in tabular form, For the purpose of illus- 


tration, a sample calculation sheet is shown in Table E-l. 
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